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Classic experiments on the distribution of ducks around separated food sources found consistency 
with the ‘ideal free’ distribution in which the local population is proportional to the local supply 
rate. Motivated by this experiment and others, we examine the analogous problem in the microbial 
world: the distribution of chemotactic bacteria around multiple nearby food sources. In contrast 
to the optimization of uptake rate that may hold at the level of a single cell in a spatially varying 
nutrient field, nutrient consumption by a population of chemotactic cells will modify the nutrient 
held, and the uptake rate will generally vary throughout the population. Through a simple model we 
study the distribution of resource uptake in the presence of chemotaxis, consumption, and diffusion 
of both bacteria and nutrients. Borrowing from the held of theoretical economics, we explore how 
the Gini index can be used as a means to quantify the inequalities of uptake. The redistributive 
effect of chemotaxis can lead to a phenomenon we term ‘chemotactic levelling’, and the inhuence of 
these results on population htness are briefly considered. 

PACS numbers: 87.18.Gh, 87.17.Jj, 87.23.-n, 89.65.-s 


INTRODUCTION 

In one of the more amusing, yet influential experiments 
on animal behavior, Harper [1] studied the distribution 
of mallards around two separated sources of standard¬ 
ized pieces of bread. After an induction period on the 
order of a minute, the average number of ducks clus¬ 
tered tightly around each station stabilized. The distri¬ 
bution he observed was simple: the number of ducks at 
each source was proportional to the flux of bread there 
(pieces/minute). This constituted the first experimental 
observation of the so-called ideal free distribution previ¬ 
ously introduced in theoretical ecology. Using the termi¬ 
nology of Fretwell and Lucas 0, ‘ideal’ means that ducks 
can identify the source where their uptake is maximized, 
and ‘free’ implies unfettered ability to access the source of 
choice. This distribution, resulting from individual ratio¬ 
nal behaviors, achieves a population-wide uniformization 
of the probability of uptake, and can be understood as a 
Nash equilibrium [3]. These works impacted not only or¬ 
nithology, but ecology HUE evolutionary biology |6j and 
the study of human behavior all areas involving 

resource acquisition in a heterogeneous environment. 

Here we take motivation from Harper’s experiment, 
and others discussed below, to examine resource acqui¬ 
sition in a heterogeneous microbial world m where 
swimming microorganisms respond to nutrient sources 
through concentration fields determined by molecular dif¬ 
fusion and microbial uptake. For the specific case of per- 
itrichously flagellated bacteria such as E. coli and B. sub- 
tilis , cells move in a run-and-tumble random walk biased 
by concentration gradients, resulting in drift of the popu¬ 
lation up these gradients nu. Because chemotaxis [l2| is 
quite different from the visually-based searching of higher 
animals, and because of the diffusive behavior of nutri¬ 
ents and the cell populations the microbial problem is 


distinct in character. This feature motivates the present 
investigation of the consequences of a collection of in¬ 
dividual chemotactic responses on the population-scale 
distribution of resources. While chemotaxis is generally 
thought to optimize uptake at the single-cell level Il3l . 
even the mere presence of translational diffusion in a pop¬ 
ulation would lead to a distribution of uptake rates. And 
the interplay of chemotaxis and consumption will mod¬ 
ify the distribution of resources and cells, with further 
potential impact on the uptake rate distribution. In this 
paper we focus on three key questions in this area: What 
is the distribution of bacteria around spatially distinct 
nutrient sources and their associated impact on the re¬ 
source field? What is the distribution of resource uptake 
rates within that population? What are the consequences 
of such distributions for cellular fitness? 

A historically important experiment on spatially- 
varying resources is Engelmann’s 1883 determination of 
the action spectrum of photosynthesis HU- Having dis¬ 
covered bacteria that are attracted to the oxygen pro¬ 
duced by photosynthesis HE he imaged the solar spec¬ 
trum onto a linear algal cell in an air-tight chamber con¬ 
taining such bacteria. They clustered around the alga 
in proportion to the local oxygen production, revealing 
with greater precision than the available techniques of 
the time the peaks of photosynthetic activity for blue 
and red wavelengths. How reliably the local bacterial 
accumulation reflects the oxygen production rate, in the 
face of both bacterial and oxygen diffusion, remains an 
open question in the spirit of the present investigation. 
Moreover, this work demonstrates how micro-domains re¬ 
leasing a limited quantity of attractive nutrients in a con¬ 
tinuous fashion can arise in the microbial world. 

Understanding bacterial organization and uptake 
around algal resources also finds an important biologi¬ 
cal context in the case of bacterial-algal symbiosis. The 
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recent discovery |16] that many algae dependent on vita¬ 
min B12 obtain it through symbiotic relationships with 
bacteria raises questions about spatio-temporal aspects 
of symbiosis: how the two species find each other and 
arrange themselves to achieve the symbiosis. 

These are examples of a more general problem of mi¬ 
croorganisms responding to the ‘patchy’ nature of nu¬ 
trients in ecosystems mi, such as the prosaically named 
‘marine snow’ [T8]. We emphasize the fundamental differ¬ 
ence between live microbial sources, such as Engelmann’s 
alga, and inert sources, such as lysis events; the former 
continuously release nutrients at low rates. They are thus 
stable in time but can be significantly impacted by bacte¬ 
rial uptake. These characteristics makes them the natu¬ 
ral microbial equivalent to the limited continuous sources 
considered by Harper in his animal experiments. 

To make concrete the interplay between production, 
consumption, diffusion, and chemotaxis we consider a 
generalized Kcller-Segel (KS) model [T9] - While the KS 
model is a well-established model that has found frequent 
application in the study of spatially extended microor¬ 
ganism populations [20], the biological setting of local¬ 
ized, low-intensity sources with a steady release of nutri¬ 
ents is little-studied, and the overarching issue of resource 
uptake distribution is essentially unexplored. In what fol¬ 
lows the steady-state distributions arising from multiple 
localized sources are analyzed to understand the conse¬ 
quences of chemotaxis on the distribution of uptake in the 
bacterial population, the total uptake being fixed. Bor¬ 
rowing from theoretical economics, we next propose that 
the Gini index [22] , a number originally used to char¬ 
acterize wealth inequality, can be used to quantify the 
individual uptake distribution. Varying the model pa¬ 
rameters and the dimensionality of space we show that 
chemotaxis can switch from redistributing the resource 
to generating greater inequalities. Finally, we explore the 
potential biological consequences of uptake redistribution 
through an example of growth at low nutrient levels. 


THE MODEL 

Consider bacteria with mean concentration bo and lo¬ 
cal concentration b(r,t) in a d-dimensional volume L d . 
Within the volume are nutrient sources with fluxes {4>i} 
of typical value </>o leading to a concentration field c(r, t). 
b and c obey KS equations m, 

— = D c S7 2 c — bf(c) , (1) 

g=D b V 2 &-V-[ X (c)fcVc] , ( 2 ) 

with D c and D & the diffusion coefficients of nutrients and 
bacteria respectively, and \(c) the chemotactic response 
coefficient. The nutrient uptake rate /(c) per bacterium 
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FIG. 1. (color online). Numerical results in one dimension. 

(a) Sketch of the setup, (b) Example of a steady-state dis¬ 
tribution of bacteria in the rescaled KS model |d| & |5]). (c) 
Corresponding distribution of nutrient at steady state. In 

(b) & (c) the parameters values are A = 3 for the domain 
size parameter and s = 1/8 for the relative strength of the 
left source. The importance of chemotaxis against bacterial 
diffusion increases with the chemotactic parameter a. 


is expected to behave at low c as /(c) ~ kc , with satu¬ 
ration at high c: /(c) ~ k ma x- When the nutrient is es¬ 
sential for life, such as oxygen for obligate aerobes, /(c) 
may vanish below some c* (24] [25] . As in Harper’s study 
and Engelmann’s experiment, the interesting regime has 
the resource limiting, with bo and L such that the to¬ 
tal nutrient flux can be consumed and is thus below the 
maximum value k max boL d . Otherwise, steady state can¬ 
not be attained and c increases indefinitely. 

For small c, we identify the length scale = 
(Dc/kbo) 1 / 2 for concentration gradients due to uptake. 
For the a run-and-tumble chemotaxis mechanism to op¬ 
erate and the continuum model to be relevant, (k should 
be large compared to the run length ( run = vt, where 
v is a typical swimming speed and r the time between 
tumbles. That is, we require the Knudsen-like number 
lrun/lk 1. For E. coli , with v ~ 20 gm/s and r ~ 1 
s, (run ~ 20 /irn, but it can be considerably longer for 
other bacteria [28]. Still at low c, with x(c) = Xo [21] . 
interesting behavior occurs when the chemotactic flux 
~ Xob 0 (/>o/Dc dominates the diffusive flux ~ -D(,6 0 /^fc, 
with (j>o/D c a typical concentration gradient. Thus, the 
Peclet-like number 


Xo4>o 

DbD c (k 


(3) 


exceeds unity. 

A one-dimensional version of Harper’s experiment 
(Fig- 0 a)) has fluxes 4>i and <f> r at the left and right 
domain boundaries, and <f> = + <p r . Nondimensional- 

izing as above we obtain in the low-c regime 

(4) 

S- 1 — =V 2 b-aV -[bVc] , (5) 
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where S = Db/D c , with the following nondimensional 
boundary conditions: 

n • Vc| 0 = — s n • (V6 — a6Vc)| 0 = 0 (6) 

n-Vc| A = l — s n ■ (V6 — a&Vc)| A = 0 , (7) 

where n is the outward unit normal to the domain. 

Two new parameters drive the evolution of the sys¬ 
tem: A = L/£k, the domain size relative to the screening 
length, and the relative strength s = </>;/$ of the left 
source. Before investigating the steady-state solutions to 
Q and (|5|, we point out that by construction the total 
uptake U over the population is equal at steady-state to 
the total nutrient flux into the chamber, 


,A 

U= dxbc 

Jo 


dxX7 2 C = n- Vein = 1 


( 8 ) 


independent of the choice of parameters and strength of 
chemotaxis. Thus overall uptake optimization is not part 
of the present study: our interest resides instead in how 
the individual behaviors result in a distribution of finite 
resource among population members, just as in Harper’s 
experiment. 

One of the natural questions to ask is whether there is 
a variational structure to the KS equations Q and [5]). 
The only related result of which we are aware concerns 
the case when the nutrient consumption rate /(c) takes 
the aforementioned high-c form of a constant. After suit¬ 
able rescaling that dynamics is 

^ = V 2 c + vb , (9) 

r)h 

<5 -1 —= V 2 & — aV • [5Vc] , (10) 

where v = — 1. As shown recently )26l . a dynamics of 
a closely related form is variational. If we introduce the 
energy functional 


£[b,c] 



— b log b 
a 


6c+ 2 |Vc| 2 


then the variational relations 


( 11 ) 


here, with /(c) ~ c also appears not to possess a vari¬ 
ational structure and thus it is not possible to conclude 
that the steady-state solutions are in any way minimizers 
or maximizers of some energy-like functional. It follows 
that the distribution of uptake rates is a nontrivial fea¬ 
ture of the underlying diffusion-consumption-chemotaxis 
dynamics. 

Figures [ljb) and[TJ c) show steady state distributions 
of b and c for a = 3, A = 3, 6 = 0.6 and s = 1/8, which 
correspond to physical values of the dimensional param¬ 
eters ( D c = 5 x 10 -6 cm 2 s -1 , Db = 3 x 10 -6 cm 2 s -1 , k = 
2 x 10 -7 cell -1 s -1 cm 3 , xo = 2.7 x 10 -3 cm 2 s -1 mM -1 , 
bo = 10 6 cells cm -3 , </o = 3.4 x 10 -6 mMcms -1 and 
L = 150 fi m pS SO]), along with the non-chemotactic 
case and the case a = 7. Bacteria accumulate on both 
sides, with more closer to the stronger source, as one 
might expect. Intriguingly, this leads to what we term 
chemotactic levelling of the nutrient: a more uniform 
concentration field than without chemotaxis. In par¬ 
ticular, we notice that the maximum uptake rate of a 
bacterium in this population, obtained closest to the 
strongest source, is decreased by chemotaxis. Diffusion 
of b and c precludes the ideal free distribution. 


QUANTIFYING INEQUALITIES WITH THE 
GINI INDEX 

Whereas studies on biological consequences of chemo¬ 
taxis usually measure the increased uptake over the 
whole or part of the population with respect to the non- 
chemotactic case (e.g. [31] [32]), as emphasized above, the 
present study is focused on the case for which the mean 
uptake rate is independent of the chemotactic behavior. 
As our interest then resides in how equally resources are 
spread among the population, comparison of the chemo¬ 
tactic results with the non-chemotactic reference distri¬ 
butions thus requires a measure of the proximity with the 
ideal free distribution. Among the many possible mea¬ 
sures of inequality I33j . we consider here the Gini index 
G £ [0,1] [223 j which for a distribution P(w) of wealth w 
in a population can be expressed as 


dc 

dt 


Sc 


r _i db 

6 a =aV ' 


bV 


SP 

Jib 


( 12 ) 


yield ([9]) and (10), but with v = +1. This case corre¬ 
sponds to the well-studied situation in which the bac¬ 
teria are sources for the chemoattractant, rather than 
sinks. Because the same term (—be) in £ yields both the 
production/consumption term in ([£]) and the chemotac¬ 
tic term in ( |10| ), the case v = — 1 appears not to have any 
variational structure of this type m- But since the con¬ 
sumption rate per bacterium is constant in this regime, 
the distribution of resource acquisition rates is trivial 
and a variational structure would provide no new infor¬ 
mation. More importantly, the case under consideration 


JJdudvP(u)P(v)\u - v\ 

2 JduuP(u) ' 1 J 

The ideal free distribution, in which every individual has 
the same wealth uq , is P = S(u — uo) and thus G = 0, 
while larger values hold for more unequal distributions. 
G can be used with any notion of ‘wealth’ jMj, such as 
biodiversity [55] , We should empasize that in using the 
Gini index to quantify uptake inequalities in the present 
study we do not imply any preferred status to G as a 
metric for resource acquisition distributions. There are 
many that could be explored, the Gini index presenting 
the advantage of being easily translated in the context of 
the continuous distribution and thus enabling analytical 
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FIG. 2. (color online). Variations of Gini index for bacterial 
uptake, with s = 0 (single source on the right side), (a) G at 
steady-state with respect to the chemotactic parameter a for 
various domain sizes A. (b) Ratio of G in the chemotactic case 
to that without, for various chemotactic strengths a , showing 
optimal domain size parameter for relative decrease of G. 


work on its variations. Using the individual nutrient up¬ 
take rate as wealth, we can transform the integrals in the 
uptake rate levels into integrals in space. This makes the 
bacterial density appear as the equivalent in space of the 
frequency distribution of uptake levels. Normalizing this 
new density function, thus making the integral for the 
total number of bacterial cells appear at the denomina¬ 
tor, we finally re-express (13) for time-dependent spatial 
distributions in a domain Q: 


G(t) = J n J n dxdyb(x,t)b(y,t)\f(c(x,t)) ~/(c(y,f))| 
2f n dxb(x,t)f(c{x,t))J n dxb(x,t) 

(14) 

With this measure of inequality, we investigate how 
the uptake distribution at steady-state depends on the 
parameters of the KS model Before proceeding 

we should make clear that the Gini index values calcu¬ 
lated within the present approach and the corresponding 
inequalities in the uptake levels are instantaneous: phys¬ 
ically, bacteria would swim along a biased random walk 
inside the steady-state nutrient distribution, thus sam¬ 
pling different concentration levels. Over a time long in 
comparison with the typical time of bacterial diffusion at 
the scale of the experimental chamber this motility would 
tend to level the integrated uptake within the popula¬ 
tion and yield lower values of G. The fundamental issue 
is then whether the time scale for this smoothing-out of 
inequalities is large or small compared to the timescale 
Ti n t for a relevant internal biological process based on 
nutrient uptake. The present approach is thus in the 
limit T int D b /L 2 <C 1, thus most relevant to large system 
sizes and short internal times. As an example, it takes 
approximately 17h for a typical run-and-tumble bacteria 
(D b = 4.10” 6 cm 2 .s _1 ) to explore the space between two 
sources separated by 5 mm, a time that is much bigger 
than the typical scale of key cellular processes such as 
division (approx. 30 minutes). 

We now go back to investigating the role of each pa¬ 


rameter: the ratio of diffusion coefficients 5 impacts tran¬ 
sient dynamics but does not modify steady state solu¬ 
tions. From numerical solutions in the phase space de¬ 
limited by s £ [0, 0.5] (from one source on the right to 
equal sources), a £ [0,15] (strength of chemotaxis) and 
A £ [0.5,15] (domain size), we obtain first the intuitive 
result that for given chemotactic and domain size pa¬ 
rameters, the more equal the sources are, the more equal 
the uptake is among the population, and the lower is 
G. More balanced sources indeed create smaller nutri¬ 
ent gradients, thus a lesser range of uptakes and weaker 
chemotaxis. We also find that G increases with the size 
parameter A (Fig. [2] (a) in the case s = 0), for larger A 
corresponds to stronger variations of the concentration 
field and higher variations of uptake. 

The question of whether chemotactic levelling of the 
nutrient concentration field can make the distribution of 
individual uptake more ideal is addressed by varying a. 
Its impact is subtle: chemotaxis levels nutrients across 
the domain, but accumulation of cells near sources im¬ 
proves uptake of some to the detriment of others. We find 
that G actually decreases with a (Fig. [2][a) for a single 
source), which reveals chemotactic levelling of the uptake 
rate among the population. In Fig. [T] G decreases from 
~ 0.3 with no chemotaxis to ~ 0.25 when chemotaxis 
(a = 3) is allowed. 

This decrease is best understood for a single source on 
one side of the domain. Because c decreases monoton- 
ically from the source, the bacterial population can be 
split into quantiles of uptake that are ordered in space. 
Chemotaxis lowers the nutrient concentration over the 
whole domain, but mostly close to the source, and it 
shifts the center of mass and mean uptake level of each 
quantile. Fig. [3jr show that bacteria with the higher up¬ 
take, which are also closer to the source, are transferred 
to lower uptake levels. For the lower uptake quantiles, 
higher levels of uptake compared to the non-chemotactic 
case are attained due to chemotaxis toward the source. 
Together, these effects bring uptake levels closer to the 
average, lowering G: the bacterial system moves closer 
to the ideal free distribution. 

The generality of this result can be established in the 
limit of weak chemotaxis (a <C 1) still with a single 
source, where a series solution G ~ Go (A) -t-ctGi(A) + • • • 
of Q and ^ yields 


G 0 (A) = 1 — 2 

a ^ = -ix 


cosh A — 1 
A sinh A 
cosh A — 1 
^ A 2 sinh A 



A cosh A A 
3 sinh(A) J 


(15) 
, ( 16 ) 


where Gi(A) < 0, so G indeed decreases with chemotaxis. 
A cumbersome analysis (not shown) for the general case 
s £ [0,0.5], that is G — Go(A, s) T ol G^(A, s) T * * *, also 
shows that Gi(A, s) < 0, thus extending this result to 
any ratio of source stengtlrs. 
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FIG. 3. (color online). Interpretation of steady-state Gini 
index variation in different dimensions d. (a) d = 1 with do¬ 
main size parameter A = 3, s = 0 (single source on right side), 
chemotactic strength parameter a = 3. (b) d = 2 with a cir¬ 
cular source of size l = 0.5 in a domain of radius L = 4.5, 
and chemotactic strength parameter a = 3. Shown are center 
of mass and mean uptake level of the sextiles of uptake dis¬ 
tribution (colored circles), overlaid on nutrient concentration, 
without and with chemotaxis. Dashed line shows concentra¬ 
tion for average uptake. Inset: G vs. chemotactic strength 
parameter a , with blue circles corresponding to displayed so¬ 
lutions. 


In the limit of strong chemotaxis (a 1), where bac¬ 
terial diffusion becomes irrelevant, we may expect to re¬ 
cover the ideal free distribution. Analytical progress on 
this steady-state problem is achieved by integrating twice 
^ to obtain 5(c), substituting into Q and then expand¬ 
ing in powers of 1/a. One obtains 


b(x) 
c{x) 


1 Aw 2 

a 2 cos 2 [x (x — j3i) /2] 

\ + — {/?2 - 2 log [cos {x(x 


(17) 

/3i)/2)]} , (18) 


where x and /3i are determined through 


x tan 


(| A ) =aS 

(—/3i + A) = a(l — s) . 


ix tan ( — 
ix 
12 


(19) 

( 20 ) 


The constant /?2 depends on x and the model parameters. 
Neglecting terms in 1/a 2 and higher order, it can be 
written as 


/32 — — 2 log a + 2 + 2 log x 

— 2 [s log s T (1 — s) log(l — s)] — x 2 A/a -|-. 

(21) 

This solution enables us to obtain an analytical expres¬ 
sion for the Gini index in the limit of strong chemotaxis: 


G~A<1 — 2s + 4 s 2 + s(l — s)log 


( 1- S ) 2 



( 22 ) 

to leading order in 1/a, for s £ [0,0.5]. This establishes 
further the generality of its decrease with a together with 
its increase with A: in ID, chemotaxis levels the uptake 
throughout the population. Moreover, in this range of 
high a, inequalities of uptake initially increase when the 
system changes from a single source to more balanced 
sources. 

As this levelling of the uptake distribution appears as 
the microbial equivalent of the more uniform uptake dis¬ 
played by ducks, it is natural to ask if, as in Harper’s 
experiments, the bacterial population reaches the ideal 
free distribution as a —A oo and splits into two localized 
sub-populations proportional to the source strengths. If 
we consider that the position Xq of minimum bacterial 
concentration separates a left population Bl associated 
to the left source and its right-hand-side equivalent H/j, 
our analytical solution for a>l directly yields 


rPi 

/ dxb(x) = Bl - As . 

Jo 


(23) 


Thus, the population associated to one source is directly 
proportional, to leading order, to the flux of this source, 
as in the ideal free distribution. Moreover, analysis of 
(181 shows that in the limit a —I oo, b{x) is localized 
in regions of width ~ 1 /y/a at both x = 0 and x = A, 
with peak values ~ a at these positions. In the limit 
of a —> oo, we thus get a localisation of the number of 
cells proportional to the source at the source: this is the 
(unphysical) limit of a microbial ideal free distribution. 

When the domain size parameter A > 1, the central 
portion of the domain has a steady-state concentration 
c ~ 0, with very small gradients. Bacteria there are 
screened from the sources, unable to feel sufficient gradi¬ 
ents to move chemotactically closer to them. The relative 
redistributive effect of chemotaxis compared to its ab¬ 
sence must then diminish with distance. Considering the 
relative Gini index G(a)/Go in the approximation (15)- 
(16), we indeed find an optimal domain size, Xq — 3.12, 
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for which the redistributive effect of chemotaxis is the 
strongest. An optimal size is also found in simulations 
beyond the linear regime in a (Fig. [2] (b)) and with in¬ 
flux from both sides, with Xq ~ 3. The decrease of this 
relative change for high values of A embodies the afore¬ 
mentioned screening, while the behavior at low A results 
from a nearly uniform concentration over the domain, 
with only weak gradients for a chemotactic response. 

Does the uptake levelling found in d = 1 hold in higher 
dimensions? To answer this, we solve Q and ([5]) for a 
single spherical source of radius l in a closed spherical do¬ 
main of radius L, both measured in units of Ik- We find 
that in d = 2 and 3 the effects of chemotaxis, for a given 
size of source and domain, are much weaker. Moreover, 
for certain parameter values, chemotaxis can actually in¬ 
crease G (inset in Fig. [3] (b)). Analysis of the quantiles 
(Fig. [ 3)3 for a two-dimensional example) shows that in 
these cases, even though bacteria closest to the source 
have a lowered uptake, a majority of the bacteria that are 
already above the average uptake in the non-cliemotactic 
case gain access to even higher uptake. Bacteria furthest 
from the source, and below the average uptake level in 
the non-chemotactic case, see their mean uptake decrease 
even further. Overall this corresponds to an increase of 
G: in higher dimensions, chemotaxis can bring the bacte¬ 
ria further away from the ideal free distribution, that is, 
it increases the inequalities among the population. The 
increase or decrease of inequalities of uptake, as revealed 
by the positive or negative change of G, may thus depend 
in detail on the system characteristics embodied in l, L 
and a. 


IMPLICATIONS FOR FITNESS 

What would be the biological consequences of chemo¬ 
tactic levelling of resources and of uptake rates? Uptake 
of nutrients governs a wide range of bacteriological pro¬ 
cesses, among which are cell growth and division. In par¬ 
ticular, the yield of biomass per unit nutrient taken up is 
an increasing function of the available nutrient concen¬ 
tration, a feature which has been suggested as selective 
for the response characteristics of chemotaxis m Here 
we provide a brief discussion of how the chemotactically- 
driven redistribution of resources throughout a popula¬ 
tion can impact on the average growth rate of the popu¬ 
lation which we consider a measure of fitness. Continuing 
the point of view taken in the introduction, we show that 
while at the single cell level in a defined resource field 
chemotaxis may increase fitness, this is not necessarily 
true at the population level. 

We compute an average growth rate /2 over the pop¬ 
ulation from the steady-state distributions of © , in a 
model in which the local growth rate y(x) is proportional 
to the local uptake rate through a yield function y(c): 


y(x) = y(c) c and 

1 f x 

fj, = — / dxy(c)cb . (24) 

* Jo 

Whereas the mean uptake is fixed by the boundary 
conditions in our problem, this average growth coefficient 
will depend on the distribution of the resource among 
the bacteria and thus will be modified by chemotaxis. In 
order to capture the increase of y with c, we consider that 
y = 0 below of threshold concentration c m i n , and adopt 
a Michaelis-Menten form above it m- 

y{c) = yoj -^-— C,n ."^ for C > Cmin (25) 

(c - Cmin) + A 

with K a saturation constant of the yield (Fig. |4^,). 

The threshold concentration c m i n can be considered as 
the limit below which all the uptake is directed toward 
maintenance costs [36j. The relative change of the aver¬ 
age growth rate due to chemotaxis ( p-chemo ~ Pnonchemo) 
is shown in Fig. B>- We observe that for lower values 
of the threshold concentration c m in the redistributive ef¬ 
fect increases the number of cells that reach the growth 
threshold: the population fitness becomes higher with 
chemotaxis than without. However, for higher values of 
the threshold concentration c m i„, the redistribution of 
uptake throughout the population leaves more cells be¬ 
low the growth threshold: we get the counter-intuitive 
result that chemotaxis effectively lowers population fit¬ 
ness with respect to the non-chemotactic case. This re¬ 
sult, which stems from the competition of the bacteria for 
the same resource, shows that in the context of contin¬ 
uous sources, an homogeneously chemotactic population 
could be selected against due to the dilution of a scarce 
resource resulting from the chemotactic behavior. 

CONCLUSIONS 

We have shown that the organization of bacteria 
around localized nutrient sources is fundamentally dif¬ 
ferent from that of higher animals due to diffusion of re¬ 
sources and feeders. Yet, there are still common charac¬ 
teristics. First, what might be termed ‘foraging’ behavior 
decreases the maximum uptake rate through competition 
for the resource; in the bacterial case this corresponds to 
a decrease of the maximum of the concentration field with 
chemotaxis. Second, foraging generates, quite evidently, 
localization and accumulation of the population closer to 
the resources. But whereas the conjunction of these two 
phenomena brings Harper’s ducks to the ideal free dis¬ 
tribution, it may fail in the microbial world: it brings 
the system closer or further from this ideal distribution 
depending on the spatial dimensionality and parameters 
capturing the strength of chemotaxis, the size of the re¬ 
sources and the distance between them. The redistribu¬ 
tion of uptake is not without consequences: when the 
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Normalised minimum concentration c m i n 

FIG. 4. (color online). Consequences of uptake levelling on 
population growth rate. Using a = 4, A = 3 and s = 0.125 
in d = 1. (a) Yield as a function of nutrient concentration, 
for Cmin = 0.2, K = 0.05 and yo = 1. (b) Change of growth 
rate as a consequence of chemotaxis for K = 0.05 and varying 
Cmin' if chemotaxis can increase growth for low yield thresh¬ 
old, the uptake leveling makes it less advantageous for higher 
threshold by diluting the resource. 


resource is scarce in comparison to the metabolic needs, 
chemotaxis effectively dilutes it and reduce the average 
population fitness. 

The issues addressed here suggest experimental stud¬ 
ies of model systems in physical ecology for which in situ 
measurements of local metabolic activity and nutrient 
concentration fields are possible. Optically-based quan¬ 
titative measures of photosynthetic activity [35], probes 
of local oxygen concentration [33] , and local mass spec¬ 
trometry m are examples of relevant techniques. Micro¬ 
bial communities in biofilms, sediments |41j . and algae 
sustaining a motile population of bacteria around them 
by releasing oxygen [52] represent interesting systems in 
which to study the distribution of uptake rates. 
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